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f"^ . Abstract. 

£SJ ' High frequency velocity oscillations were observed in the spectral lines Fe I 543.45nm 

and 543.29nm, using 2D spectroscopy with a Fabry- Perot and speckle reconstruction, at 
the VTT in Tenerife. We investigate the radial component of waves with frequencies in the 
h^H ' range 8 - 22mHz in the internetwork, network and a pore. We find that the occurrence of 

, waves do not show any preference on location and are equally distributed over down-flows 

■ and up-fiows, regardless of the activity of the observed area in the line of Fe I 543.45nm. 

The waves observed in the lower formed line of Fe I 543.29nm seem to appear preferentially 

1 i , over down-flows. 

6 '■ 
to: 

1. Introduction 

C3 . 

, According to Kalkofen (2001) the problem of chromospheric heating is a 

question of energy supply for the radiative emission. It is thought that 
high frequency acoustic oscillations may provide this energy supply. The 
\ generation of acoustic oscillations can be described by the 'Lighthill mech- 

I ■ anism' (Lighthill, 1951; Proudman, 1952). Deubner (1983) argues that high 

. frequency bursts are either generated by rising granules, or propagate more 

| or less uniformly from deeper layers into the convection zone. 

The frequencies of acoustic oscillations depend on various fluid flow pa- 
rameters (Lighthill, 1951). The most energetic oscillations should be gen- 
.^J ' erated in those regions where the convective velocities are largest (Stein, 

1967). 

Kalkofen (1990) suggests that the location of acoustic oscillations depend 
on their frequency. In areas of strong magnetic fields at the cell boundary, 
Kalkofen calculates that heating is done by oscillations with frequencies of 
1 — 3mHz. The bright points are heated by oscillations with periods around 
6mHz; while areas free of magnetic fields will be heated by oscillations of 
higher frequencies. 

Observations (Espagnet et al, 1996) show that acoustic events occur 
preferentially in dark inter granular lanes, i.e. corresponding to down flows 
of plasma and conclude that the excitation of solar oscillations is associated 
with rapid cooling occurring in the upper convection layer. 
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Dominguez (2004), Andjic and Wiehr (2006), Andic and Vockic (2007) 
and Trujillo Bueno et aZ.(2004) have shown evidence for magnetic activity in 
the quiet sun. Domingues (2004) found that the magnetic flux in the quiet 
internetwork is smeared by insufficient resolution and therefore underesti- 
mated. He also found that there are areas in the quiet internetwork which 
have magnetic fields comparable to those from the pores. Trujillo Bueno et 
al. (2004) claim that we can only see 1% of the magnetic field present in 
quiet internetwork regions. Andjic and Wiehr (2006) and Andic and Vockic 
(2007) show evidence for some magnetic activity in the quiet internetwork 
regions used in this work. 

2. Observations 

The instrument used in this work is the German Vacuum Tower telescope 
(VTT) on the Canary islands, with the Fabry-Perot spectrometer (Bendlin, 
Volkmer, and Kneer, 1992). The instrumental setup has been described in 
detail by Koschinsky et al. (2001). We have not used a polarimeter in this 
study. The observations were performed using the spectral lines of Fe I at 
543.45nm and 543.29nm in June 2004. The data sets used in this work are 
listed in Table I. 



Table I. The data sets used in this work, dates, used lines and object of observation. Position 
[0, 0] corresponds to the Sun center. 



Mark 


Lines [nm] 


Area 


Coordinates 


Images 


Expo- 


Cad- 


Dur- 












sure 


ence 


ation 












[ms] 


[s] 


[min] 


DS1 


543.45 & 543.29 


Quiet Sun 


[0,0] 


108 


30 


28.4 


52.7 


DS2 


543.45 & 543.29 


Bright points 


[0,0] 


108 


30 


28.4 


33.61 


DS3 


543.45 


Pore 


[0,0] 


108 


30 


28.3 


23.19 


DS4 


543.45 


Quiet Sun 


[96.7, 90.7] 


91 


20 


22.7 


52.9 



3. Data Reduction and Analysis Methods 

The procedures used for data reduction in this work are the same as those 
used in Andic (2007) and Andjic (2006). 

The data reduction for the broad-band images includes the following 
procedures: subtraction of darks, flat fielding and speckle reconstruction. 
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The reconstruction is performed with software developed by P. Sutterlin 
following the method of de Boer (1992, 1994). The program used for the 
reconstruction of narrow-band images was developed by Janssen (2003) and 
is based on the method developed by Keller et al. (1992). The velocity maps 
were calculated for the line cores from the line bisectors using the technique 
described in detail by Andic (2007). Line profiles were calculated using 3D 
radiative- hydro dynamical model of Asplund et al. (2000). This method is 
described in detail in the work of Shchukina and Trujillo Bueno (2001). The 
formation height in LTE is 308.2km for Fe I 543.29nm and 660.2km for Fe 
I 543.45nm. In NLTE the corresponding heights are 258.5km and 588.7km 
for Fe I 543.29nm and Fe I 543.45nm respectively (N. Shchukina, private 
correspondence) . 

For the wavelet analysis used in this work we chose the Morlet wavelet: 

fo(t) = 7 r- : *e i " ot e-%, (1) 

where lvq is the non-dimensional frequency and t the non-dimensional time 
parameter. The wavelet analysis code used is based on the work by Torrence 
and Compo (1998). The associated Fourier period, P, is 1.03s for ujq = 6, 
where s is the wavelet scale (see Table 1 in Torrence and Compo (1998)). 
One-dimensional wavelet analysis is done for each spatial coordinate. At the 
beginning and end of the wavelet transforms there are regions where spurious 
power may arise as a result of the finite extent of the time series. These 
regions are referred to as the cone of influence (COI), having a temporal 
extent equal to the e-folding time (td) of the wavelet function. In our case it 
is td = V%s = • This time scale is the response of the wavelet function 

to noise spikes and is used in our detection criteria by requiring that reliable 
oscillations have a duration greater than td outside the COI. This imposed 
a maximum period of 912s (or « lmHz) above which any detected periods 
were disregarded. The automated wavelet analysis carried out in this work 
has previously been presented in detail by Bloomfield et al. (2006). In the 
current study only the power outside the COI is retained. 

To determine whether registered oscillations are located above down-flows 
or up-flows, a comparison with velocity maps is performed. In this analysis 
the power has been normalized and only power in excess of 15% of the 
normalised value is considered (Figure 1.). 

A comparison of the spatial distribution of high frequency oscillations 
with white-light structures was performed. 

High frequency oscillations power maps are compared pixel by pixel with 
the corresponding white-light image. In Figure 2. a scatter plot for one frame 
in the time sequence is shown. The temporal difference between the white- 
light image and the power map is calculated according to the assumption 
that the high frequency oscillations travel with the speed of sound. 
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Figure 1. Velocity map from one of the VTT data sets, where the field of view shows only 
those velocities which appear at the same location as the high frequency oscillations, from 
the frequency interval 22mHz to 8mHz, with power above 15% of the maximum level. The 
remainder is set to a zero value. 



This assumption is oversimplifying the true picture of the lower solar 
atmosphere and the oscillations within. Up to the present, theoretical models 
of stellar atmospheres have been based on several simplifying assumptions, 
such as the restriction to a 1-dimensional (ID) stratification in hydrostatic 
equilibrium with only a rudimentary parametrization of convection through. 
Based on such models (e.g. Howleger and Miiller 1974) it is assumed that in 
the solar atmosphere the acoustic transit time is approximately 5 minutes 
(for a height of 2000km and a sound speed of 7kms -1 )). For the propagation 
of acoustic waves in the solar atmosphere over several scale heights, their 
frequency has to be above the cut-off frequency Equation (2). 



^ac = — , (2) 

2c s y J 

where g is the gravitational acceleration, 7 = c p /c v the adiabatic coefficient, 
and c s = \/Hjg with H as density scale height, given as: 



dr 

(Stix (2002), chapter 5.2.4). Therefore in a ID model the cut-off frequency 
varies only with atmospheric height. Thus, at a certain height in the at- 
mosphere a reflection layer exists where the values of g and c s yield the 
appropriate cut-off frequency. For a given frequency, a wave can be oscil- 
lating at one height and be evanescent at another. This situation can cause 
standing waves for almost the whole range of acoustic frequencies. Fleck 
(1989) explain that there is a possibility for standing waves, originating 
from the total reflection of upward propagating waves at the chromosphere 
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- corona transition region. This discovery of standing patterns was confirmed 
by Espagnet et al. (1996). 

Virtually all information about the Sun is deduced from its emergent 
spectrum. Directly observable manifestations of Solar convection (granula- 
tion) thereby influences the emergent spectrum. Also most of our present 
empirical knowledge of solar surface magnetism stems from the analysis 
of the Zeeman effect in spectral lines. Routine magnetograms give the im- 
pression that cell interiors are non-magnetic, while some recent research 
showed that's not case (Domingues, 2004; Andjic and Wiehr, 2006; Andic 
and Vockic, 2007; Trujillo Bueno et al, 2004). 

In light of this information, the assumption that the observed oscillations 
are purely acoustic is oversimplifying and does not give the full picture 
of the processes. Never the less, since current theories for the heating of 
the chromosphere favor acoustic oscillations, part of the current analysis is 
conducted with this assumption. 

For this part of the analysis the intergranular lanes (darker areas) were 
defined as areas where the intensity of the corresponding pixel is less than 
the modal intensity of the whole field of view. 

To obtain the percentage of power which is located above dark structures, 
the power was normalized, and only that power which exceeds 15% of the 
maximum value is considered. This procedure is performed throughut the 
time sequence yielding the variation of positions with time. An example is 
shown in Figure 3. 



4. Results 

Above described analysis is performed for all data sets. Results are presented 
in Table II. The vb present the upward velocities and vr downward veloci- 
ties, while T,v presents the sum of all velocities used for this calculation. The 
formation heights for line cores are calculated for LTE and NLTE (i^LTE and 
^nlte) respectively. The percentage is calculated over the whole data set. 

The results show in Table II show that the percentages of upward and 
downward velocities are similar. This is also visible with results for T,v where 
values are less than 1% of the maximum velocity used in calculations. To 
see the time variation of this percentage averaging over the field of view is 
performed. An inspection of Figure 4 shows that there is no noticeable trend 
in the observed velocities. 

Analysis of velocities is done at each frequency. Figure 5 shows that with 
frequency there is no significant difference in the location of power , since 
there is only 0.3% variation between positions. 

Analysis of locations using intensity variations of images taken in the 
continuum of the lines showed similar results for the spectral line 543.45nm. 
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Distribution of power at core of the line Fe I 543.45nnn 




Figure 2. Scatter plot from the DS2data set. High frequency oscillations with power 
above 15% of the maximum level is presented vs. the corresponding pixel intensity in the 
granulation image. The power was integrated over the frequency range 7mHz to 17rnHz 
for this data set. The results from the spectral line 543.45nm are shown. 




Figure 3. Variation of the percentage of observed oscillations located above intergranular 
lanes from the data sets DS1 and DS2 for the Fe I 543.45nm. The solid line represents the 
percentage for the frequency range 7mHz to 17mHz from the data set DS1. The dashed 
line represents percentage for data set DS2 and the doted line marks 50% level. 



polozaj3.tex; 1/02/2008; 6:12; p. 6 



Energy of high frequency waves in the low solar chromosphere 



7 



Table II. The percentage of blue and red velocities appearing at the observed 
locations of high frequency oscillations. 



mark 


line [nm] 


VB [%] 


VR [%} 


Eu[km/s] 


Hlte 
[km] 


Hnlte 
[km] 


power 

coverage 

[%] 


DS1 


543.45 


32.22 


31.32 


-0.0086 


660.2 


588.7 


99.9 


DS1 


543.29 


11.41 


16.27 


-0.0508 


308.2 


258.6 


70.7 


DS2 


543.45 


30.24 


29.47 


0.0022 


660.2 


588.7 


100 


DS2 


543.29 


8.54 


10.75 


-0.0204 


308.2 


258.6 


81 


DS3 


543.45 


7.18 


6.6 


-0.0076 


660.2 


588.7 


72.01 


DS4 


543.45 


5.37 


5.94 


0.0031 


660.2 


588.7 


99 



variation of the dopier velocities 
0.6 I I I I 1 1 




-0.6 I I I I I 

10 20 30 40 50 60 

timefmin] 

Figure 4- Time variation of the percentage of observed velocities located above spatial 
positions for the observed oscillations. Graph is divided into two parts. The top part 
shows upward velocities and in bottom part the downward velocities are shown. 

There is no preference in location of the observed oscillations for the 
spectral line 543.45nm, as it can be seen in Figures 2 and 3. Figure 3 shows 
that the percentage of oscillations above darker areas is around 55%, which 
is not to be taken as certain indicator of a preferable position. Results for 
the spectral line 543.29nm show a tendency to appear above darker areas in 
the continuum image as it can be seen is Figures 6 and 7. In Figure 7 the 
percentage of oscillations above darker areas is around 70%. Two sudden 
drops in percentage for the data set DS2 are caused by a decrease in quality 
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-0.6 I i i i ~l ~ t i i | i i i | i i , | i i i | i i i 

6 8 10 12 14 16 18 

frequency[mHz] 

Figure 5. Frequency variation of the percentage of observed velocities located above spa- 
tial positions of the observed oscillations. The solid line represents the difference between 
upward and downward velocities (vb — vr) as a percentage for the core of the 543.45nm 
spectral line and dashed line for the 543.29nm. 



of the speckle reconstructed images due to a significant degradation of the 
seeing condition during observations. 

The percentage of the field of view which contains oscillations is given 
in Table II. This coverage is calculated by summing the positions of all 
oscillations with power above 15% of the maximum power over the complete 
time sequence. Inspection of the table shows that oscillations observed in the 
543.45nm line, which is formed higher in the atmosphere, cover more of the 
field of view than the 543.29nm line. This percentage is lower for DS3. This 
is because no oscillations were registered over the pore. 



5. Discussion and Conclusions 

The observation used in this work were obtained in quiet solar regions, active 
regions with G-band structures and an active region containing a pore. Our 
findings do not show a dependency on the selected solar locations. 

In the data sets from the VTT, the higher formed spectral line 543.45nm 
show no preference in oscillation site location (Figures. 2 and 3 and Table 
II.), while for the lower formed line, 543.29nm, there is a slight preference 
towards down flow regions (Figures. 6 and 7 and Table II). 

The influence of atmospheric seeing is mainly visible during speckle re- 
construction. The work by von der Liihe (1984, 1993, and 1994) estimates an 
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Figure 6. Scatter plot from the DS2 data set for 543.29nm. High frequency oscillations 
with power above 15% of the maximum level is presented vs. the corresponding pixel 
intensity in the granulation image. The power was integrated over the frequency range 
7mHz to 17mHz for this data set. 



Variation of the percentage of the power above the intergranular lanes 

80 I "' 1 




40 h 1 

10 20 30 40 50 60 

Time [min] 

Figure 7. Time variation of the percentage of observed oscillations located above inter- 
granular lanes from DS1 and DS2 for 543.29nm. The solid line represents the percentage 
for the frequency range 7mHz to 17mHz from the data set DS1. The dashed line represents 
the percentage for the data set DS2 and the dotted line marks 50% level. 
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error of 10% in intensity introduced by speckle reconstruction. Noise filtering 
was done following the recommendations of de Boer (1996), so that error 
levels remain at 10%. Seeing influences are most visible with differing data 
cadences. In the case of good seeing, the cadence was lower than in the case 
of very good and excellent seeing conditions. Therefore the amount of noise 
remaining after speckle reconstruction tended to vary only inside predicted 
error levels (von der Liihe, 1984; von der Liihe, 1993; von der Liihe,1994; de 
Boer, 1996). A full description of the errors is presented by Andic (2007) 

5.1. Fe i 543.45NM 

Observations (Espagnet et. al 1996) show that acoustic events occur prefer- 
entially in dark intergranular lanes. Wunnenberg et. al (2002) claim that high 
frequency oscillations occur preferentially above down-flows. Observations 
with the same spectral line as in the work of Wunnenberg et. al (2002) 
shows no evidence for this preference. One of the main differences between 
the work of Wunnenberg et. al (2002) and the present study has to do 
with the velocity response functions of the Fe I 543.45nm line. In this work, 
velocity response functions were not used. 

Figures 2 and 3 show that oscillations registered in this line do not show 
preference in location. 

Kalkofen (1990) suggests that the location of acoustic oscillations should 
depend on their frequency. However, our findings do not agree with this pre- 
diction. The difference in the location of oscillations of various frequencies is 
too small to be taken as a serious indicator of their different positions. Figure 
5 shows that differences in position with frequency do not exceed 0.3%. Also 
for this line those differences are mainly positive indicating that upward 
velocities are dominant. Only in the frequency range 9mHz to 11.5mHz are 
there indications for dominance towards donward velocities. Nevertheless, 
since the difference is less than 1% this can not be considered as a significant 
result. 

Also, all observed oscillations tend to follow similar patterns, which con- 
tradict the statement by Espagnet et. al (1996) that most energetic oscilla- 
tions are associated with down-flows in dark areas which are well separated 
from each other in time and space. Considering the fact that Espagnet et. al 
(1996) were working with oscillations of much lower frequencies, this could 
explain the difference. 

The formation height of this spectral line places the core of the line in 
the vicinity of the temperature minimum. This opens consideration that 
some of the granulation motions could be registered as oscillations in this 
work. Espagnet et. al (1995) show that intensity fluctuations associated 
with granules disappear in a very short distance, below 100km above the 
continuum level. Yet, velocity fluctuations associated with granulation cross 
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the whole thickness of the photosphere. This is typical for larger granules 
(2" — 3") above which velocities spread throughout the whole photosphere. 
Therefore it is possible that some of the oscillations detected here from 
the velocity maps could be originating from granulation. To remove such 
velocities, without imposing the limit on the registered frequency range, 
those oscillations with height depended amplitude decay are eliminated. The 
maximum velocities from averaged granules on this height were noted to 
be below lkm/h, while the maximum observed velocities in this work were 
around 6km/h, therefore all velocity oscillations with amplitudes below 15% 
were not analysed in this work. This method does not exclude velocities 
associated to very large granules where the decay is slower. 

5.2. Fe i 543.29NM 

The slight preference for oscillatory appearance above downwards velocities 
is found in observations with the this line. Comparing to the intensity of the 
continuum image shows substantial preference toward darker areas in the 
granulation pattern. Figures 6 and 7 show that oscillations registered in this 
line do show preference in their location. 

Figure 5 shows that difference in position with frequency does not exceed 
0.3%. For this line those differences are mainly negative indicating that 
downward velocities are dominant. But since the difference is less than 1%, 
this also can not be taken as a significant result. 

The formation height of this spectral line places the core of the line in 
the middle of the photosphere. This makes the results even more sensitive to 
granulation overshooting. Again, the main concern are velocities which on 
this level can be, on average, around lkm/h. Since maximum velocities noted 
for this height are also around 6km/h, the method mentioned above cannot 
successfully remove all velocities which could be associated with granulation. 
Therefore the error in the results for this line is larger, but no more than 
20%. 
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